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When two-photon interactions are induced by down-converted light with a bandwidth that exceeds 
the pump bandwidth, they can obtain a behavior that is pulse- like temporally, yet spectrally narrow. 
At low photon fluxes this behavior reflects the time and energy entanglement between the down- 
converted photons. However, two-photon interactions such as two-photon absorption (TPA) and 
sum-frequency generation (SFG) can exhibit such a behavior even at high power levels, as long as 
the final state (i.e. the atomic level in TPA, or the generated light in SFG) is narrowband enough. 
■ This behavior does not depend on the squeezing properties of the light, is insensitive to linear losses, 

f"*) ' and has potential applications. In this paper we describe analytically this behavior for travelling- 

wave down-conversion with continuous or pulsed pumping, both for high- and low-power regimes. 
£SJ 1 For this we derive a quantum-mechanical expression for the down-converted amplitude generated by 

an arbitrary pump, and formulate operators that represent various two-photon interactions induced 
by broadband light. This model is in excellent agreement with experimental results of TPA and 
. SFG with high power down-converted light and with entangled photons [Dayan et al, Phys. Rev. 

Lett. 93, 023005 (2004), 94, 043602, (2005), Pe'er et al., Phys. Rev. Lett. 94, 073601 (2005)]. 
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I. INTRODUCTION 



Parametrically down-converted light [l|, 0, H, B| exhibits correlations between the down-converted modes, and 
has been a primary source in experimental quantum-optics [j| @, 0, H, S E3, EJ- At low photon-fluxes, these 
correlations are inherently nonclassical, and are exhibited in the generation of entangled photon-pairs [HI, [l3|, [3, [l5| . 
Two-photon interactions induced by entangled photons are expected to exhibit nonclassical features, in particular a 
linear dependence on the intensity [H, [l?], [H, Il9l l20l. I2H [Hj], as was observed with two-photon absorption (TPA) 
and with sum- frequency generation (SFG) [24| . At higher powers, the correlations between the down-converted 
modes are better expressed as quadrature correlations. In this regime, the nonclassical nature of the correlations 
lies in the fact that their pre cision can exceed the vacuum shot-noise level, a phenomenon that is named squeezing 
[T3L I25I l2ol I27I [H, [29|, l30l l3ll [3^, and is expected to effect two-photon interactions even at photon- fluxes which 
exceed the single photon regime [lfj [22|, (3J, [35|, [361, [37]]. However, while squeezing (i.e. the high precision of the 
correlations) is easily destroyed (for example, by linear losses), the remaining correlations still have a dramatic 
effect on two-photon interactions. Specifically, when the down-converted bandwidth is significantly larger than 
' the pump bandwidth, the down-converted light can induce two-photon interactions with the same efficiency and 
sharp tem por al behavior as ultrashort pulses, while exhibiting high spectral resolution as that of the narrowband 
pump (38l. l39l. l40l l4lj . Although at low photon fluxes these properties are a manifestation of the time and energy 
• i-H , entanglement between the down-converted photons, at high powers exactl y th e same properties are exhibited if the 
final state of the induced two-photon interaction is narrowband enough |40j . This effect occurs since the narrow 
bandwidth of the final state "post-selects" only the contribution of photon pairs that are complementary in energy. 
Thus, in the high-power regime, these seemingly nonclassical properties are completely described within the classical 
framework (in fact, they are equivalent to the correlations that exist between the transmitted signal and its key in 
spread-spectrum communication systems (42|), and do not depend on the squeezing degree of the down-converted 
light. This equivalence of broadband down-converted light to coherent ultrashort pulses also implies an ability 
to coherently control and shape the induced two-photon interactions with pulse-shaping techniques, although the 
down-converted light is neither coherent, nor pulsed. This effect was demonstrated both in high power [39l |40| and 
with broadband entangled photons 41], in which case it can be viewed as shaping of the second-order correlation 
function of the entangled photons g 2 (r). The fact that these unique properties can be exhibited at high powers, with 
no dependence on nonclassical features of the light (such as entanglement or squeezing) , makes this phenomena both 
interesting and applicable [12, [H| . 

In this paper we describe analytically this unique behavior of two-photon interactions with broadband down- 
converted light. First, in section UTJ we analytically solve the equations of motion for broadband parametric 
down-conversion with arbitrary pump (continuous or pulsed) , obtaining a non-perturbative solution that is valid both 
for low-power and for high-power down-converted light. The derivation takes into account the specific spectrum of 
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the pump, assuming only that it is significantly narrower than the down-converted spectrum. The solution, namely 
the annihilation operators for the down-converted fields, is represented as a function of the down-converted spectrum, 
a property that is typically easy to calculate, estimate or measure. 

In section UTTl we formulate quantum- mechanical operators that can represent multiphoton phenomena, specifically 
TPA, SFG and coincidence events, induced by any broadband light. By evaluating the expectation value of these 
operators, using the annihilation operators derived in section [IIJ we obtain in section ITVl a generic analytic expression 
for two-photon interactions induced by broadband down-converted light, as a function of the pump spectrum and the 
down-converted spectrum. In section [V] we derive specific expressions for TPA, SFG and coincidence events induced 
by broadband down-converted light, and analyze their temporal and spectral properties. In section IVJ we give a brief 
summary of our results and add some concluding remarks. 



II. SOLVING THE EQUATIONS OF MOTION FOR BROADBAND DOWN-CONVERSION WITH AN 

ARBITRARY PUMP 

In our derivations we have chosen to use a continuous-variables version of the formalism suggested by Huttner et 
al. fjij . This formalism represents the quantum fields in terms of space-dependent spectral mode operators a(oj, z), 
instead of time-dependent momentum mode operators a(k,t). The advantage is that unlike momentum modes, 
temporal modes are unchanged by a dielectric medium, reflecting the physical fact that while the energy density of 
the fields depends on the medium, the energy flux does not. Accordingly, a temporal periodicity, instead of a spatial 
periodicity is assumed in the quantization process. We have performed the transition to continuous variables (taking 
the temporal periodicity to infinity) following the guidelines of the same procedure in momentum and space pB| . 
Specifically, assuming a single polarization and a single spatial mode we may write: 

E(t,z) = E + {t,z) + E~{t,z) , 

p OO I fc 

E+(t,z) = tj^ du;^/^^a(u,z)e-^ t -^ c \ E-(t,z) = {E+(t,z)Y (I) 

where S is the beam area, e is the permittivity, and a){ui, z), a(ui, z) are the slowly varying complex amplitudes of 
the creation and annihilation operators of the electromagnetic field [46}] . A major advantage of this formalism is the 
relative convenience at which we can define a momentum operator G(z) for the electromagnetic field in a dispersive 
medium, and use it as the generator for space propagation. The equations (and hence the solutions) derived in the 
following resemble those attained by using the classical Maxwell equations, or by using the Hamiltonian for time 
propagation and replacing the spatial coordinate z by tc (see, for example, [47j): however the derivation presented 
here is a multimode, continuum-frequency one, enabling us to take into account the specific spectrum of the pump 
and the fact that each signal mode is coupled by the pump to multiple idler modes, and vice versa. 



The momentum operator related to the second-order nonlinear polarization is of the form [4J, [48|, |49|, |5C 

/■OO /.OO 



/ hlJpUj(ujp — uj) . . 

P(u) p , uj) = x w, Up) J — 3^ — - — r — -r^—, r , (2) 

y 1&-Ke c 6 bn p {uj p )n 8 {uj)n l {uj p - uj) 

where the subscripts s, i,p denote the signal, idler and pump modes, respectively, n St i^ p {uj s ^ tP ) denote the corresponding 
indices of refraction, and with: 

l±.k{u) p ,u)) = k p (ui p ) — k s (ui) — ki(u> p — uj) . (3) 

The primary assumption throughout our calculations is that the down-converted spectrum A DC is considerably 
broader than the pump bandwidth 8 p : 

A DC > S p , (4) 

where we take A DC to denote the bandwidth of the signal (or the idler) field (note that the signal and the idler have 
the same bandwidth). Having assumed a narrowband pump, we can safely neglect the dependence of x (which is 
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typically real) and Afc on w p , since typically the crystal's nonlinear properties vary only at much larger scales of 

frequencies, and write: /3(u> p ,uj) = /3(w), and Afc(w p ,o->) = Ak(uj). We thus write the equation of motion for a s (u>, z) 
for undepleted, strong pump, replacing a p (u> p , z) with the spectral amplitude of the classical pump field A p (uj p ): 

±—L = --[a s (uj,z),G nt (z)] 

= -0(uj)e- iAk ^> / duu p k p {uu p )a\{uu p -uj,z) (5) 
Jo 

The same can be formulated for the idler, which leads to: 



d 2 a e (uj,z) 



, = |/3(w)| / dv duj p A p (u> p )A*(Lj p + v)a s (uj + v,z) 

OZ J-oo JO 

- Ak(uj) (3(uj)e- lAk ^ z f™ dw p A p (u p )a\(u p -u,z), (6) 



where v — top — uip. For a stationary light field we can write [51] ] : 

du) p Ap(wp) A*(uj p + v) = 2nl p 5(v) (7) 

where I p = (l p (t)) = (\A p (t)\") is the mean power (in units of photon-flux) of the pump field, and so: 

/>oo 

dv / duj p A p {uj p ) A*(cj p + v) cl s (uj + v, z) = 2nl p a s (u>, z) . (8) 
Jo 

A similar result can be obtained, for non-stationary, pulsed light with a final duration of r p , by approximating: 

f°° ^ f 7T 

/ dcjp A p (u)p) A*(uj p + v) w / I p (t) dt when \v\ < — 
Jo Jt p t p 

w otherwise . (9) 

This approximation reflects the fact that the spectral amplitude of a finite signal can be considered roughly constant 
within spectral slices that are smaller than 2tc/t, where r is the signal's duration. Accordingly, since (under the 
condition of Eq. ^) pulsed down-converted light has a duration which is always equal or shorter than the pump 
pulse, its spectral components may be regarded as constant over spectral slices that are smaller than 2tt/t p . For this 
reason it is safe to approximate: 
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dv / dujp A p (uj p ) A*(lu p + v) a s (oj + v, z) w / I p (t) dt (3(u>) / dv d s (uj + v, z) 

OO J J T B J — 

' 9 f Tp 

Ip(t) dt (3{uj) I dva s (u),z) 



= 2TrIp[3(uj) 2 a s (u,z) (10) 

obtaining the same result as with a stationary pump, except for the fact that the averaging I p — (l p (t)) is performed 
over the duration of the pump pulse. Thus, both for stationary and non-stationary pumps, Eq. © becomes: 

p*2 « ( \ />oo 

a °Y Z> = 2nl p /3(u;) 2 a s (cj, z) - Afc(w) f3(cj) e - lAk ^ z / du> p A p {u>p) a\{w p -u,z). (11) 
oz j 

The solution of Eqs. (O and (TTT]) is: 

a s {uj,z)e lAk ^ z ' 2 = ( cosh (kHz) + sinh (k(u)z) ) a.{u,Q) 

V 2k{uj) J 

Z/3(uj) /*oo 

— — -. — r- sinh (k(w)z) / duip A p (ujp) d\(u! p — uj,0) . (12) 
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with: 



k(u) = ^2ttI p (3(cu) 2 - Afc(w) 2 /4 . (13) 
The average photon flux spectral density for the signal and the idler fields is therefore [45| : 

1 f°° 

n s ^(uj,z) = — duj'(vac\al i (uj,z)a s ,i(uj / ,z)\vac) 
^ Jo 

MM 2 -:-,2 



k(w 2 ) 



sinh 2 (k(w)z) . (14) 



Since the down-converted spectrum can readily be calculated or measured for any specific down-conversion appa- 
ratus, we find it convenient to present our following calculations using n s i (w) as given parameters, thus avoiding the 
issue of evaluating and the other elements of /3(w) that determine the down-converted spectrum, and focusing on 
the behavior of two-photon interactions induced by such a light. Assuming that the down-conversion process occurred 
along a distance L, we denote: 

/ 2 \ 

Tl S,i{ U) ) = T Q( \5 n *A U > L ) = Smh2 ( K M Z ) ' ( 15 ) 

thus obtaining the following simple expression for a s ^(tu,L): 

a s>i (w,L)e iAk ^ L / 2 = Ul + n s ,(c) + ^gl^/^) \ a s >,0) 

- i ; Jn Sl i(oj) / dujp A p (ujp) a\ s (uj p - oj,Q) . (16) 



For good phase-matching conditions (i.e Afc(w) <C Wlp/3(uj)), we get k(w) — > ^J2-kI p (3(uo) and therefore 
Ak(uj) / k(ui) <C 1, in which case Eq. (fll)|) can be further simplified to: 



a aii (w, L) e lAk ^ z ' 2 = sJl+n Sti {uj) &. ti (w, 0) - — ^= y/n^(cj) / dco P A p (w p ) at s (u; p - u, 0) . (17) 



with: 

u s ,i(w) = 27rn Sii (a;,i) , (18) 

III. DERIVING OPERATORS FOR WEAK TWO-PHOTON INTERACTIONS INDUCED BY 

BROADBAND FIELDS 

In this section we derive expressions for weak (perturbative) two-photon interactions induced by broadband fields. 
Limiting ourselves to low efficiencies of interaction, we neglect the depletion of the in-coming fields a\^(j-u,z) that 
induce the interaction. Therefore we suppress in the following their dependence on z, denoting ai 2{u, z) — ai 2(w, 0) = 
ffli,2(w). We begin with specific expressions for SFG and TPA, and then obtain a generic expression which will be 
used in the following sections. 

A. SFG 

For SFG, we use the nonlinear momentum operator of Eq. @, replacing a p with the creation operator of the SFG 
mode al FG : 

da SFG (Q, z) i r /W/_yi 
q- = -- |_a SFG (S2, z), G (z)\ 

pOO 

= -i dw/3(w,ft)ai(w)a 2 (ft-w)e^ Afc( "'° )z , (19) 
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with 



Afc(fi,w) = k SFG (fl) - ki(w) - fc 2 (17 - u) . (20) 
This leads to the following approximation for a SFG (17, z = L), where L is the overall length of the nonlinear medium: 

r°° / e -iAk(u,fl)L _ ]\ 

a SFG (17, L) a SFG (17,0) + / I — — — — ) f3(aj,Q)ai(ui)a 2 (Of - u) 

Jo V A/c(u,17) J 

= a SPG (Q, 0) - 1L / due- lAk(u > n)L/2 smc[Ak(uj,n)L/2]P(uj,n)a 1 (uj)a 2 (n-uj) , (21) 



where sinc(a;) = sin(x)/x. Since we are interested only in the nonlinearly generated amplitude, we shall ignore the 
first term, d SFG (f2, 0), in Eq. (|2"Tj) . The overall SFG photon flux through the plane z = L for a given initial light state 
lib) is therefore: 



N SPG (i,L) = — / (Kl / dO / (^|a s t FG (fi,L)e ia *a SFG (r!',L) £ 
2tt Jo Jo 

= ml G (t,L)r, SF c(t,L)\ib} (22) 



where fj SFG (t, L) is the photon flux amplitude operator: 
1 f°° 

r% FG (t 7 L) = -= / dna S¥G {^L)e- int 



2tt 
ii 



dVLe- int I dLue- lAk{u > n)L/2 smc[Ak(uj,n)L/2]P(uj,n)a 1 (uj)a2{n-uj) 
Jo 

I />OC 

dfie" im / dw$(w,17)ai(w)a 2 (17-w) , (23) 



V2tt 

taking <&(w, 17) to include the coupling coefficient (5{ui, 17) and phase matching terms: 

17) ee e -«Afc(^,o)L/2 sinc [ Afe ( W) n)L/ 2 ] /3( w , 17) . (24) 

Using Taylor expansion about the center frequency of the SFG spectrum 17o, Ak(u>, 17) can be separated into two 
terms, where one depends on 17 — 17o and the other onu- (lu): 

+ ( 5^_a^ )( „_ H) + ( a^)_a^) )( _ M)2 

= Afc(17, (w))+Afc(w,f2 ) , (25) 



where the first two terms represents the difference between the group-velocities and between the group-velocity 
dispersions of the pump and the idler, and the second represents the difference between the group-velocities and 
between the group-velocity dispersions of the signal and the idler. In type-I phase-matching this implies a linear 
dependence on lu p versus a much weaker quadratic dependence on w, since the group- velocities of the signal and the 
idler arc identical: 



a,/ /<%sfg(17) <% 2 (17 — u))\ _ ^ , /d k 2 (n — lu) d ki(uj)\ , . . .„„. 

a*( W) 17) « {-^J. L_^) (17 - no) + ( b - -&W (w - • (26) 

Since the sine function in Eqs. (|2ip - (l24]) results from integration over the exponent e lAfe , and since /3(w, 17) depends 
very weakly on 17, for good phase-matching conditions (i.e. small Afc), the approximation Afc(u,17) w Afc(17, (u>)) + 
Ak(cj, 17q) enables us to represent the dependence of Q(ui, 17) as: 



$(w, 17) « 5s FG (i7 - n ) /sp G (w, n ) , 



(27) 
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with 

g SFG (Q - Oo) = e -^(n-n )L/2 sinc [ Afc ( fi _ n )L/2] 

f SFG (uj, n ) = e -^fe(^,n )L/2 sinc [ Afe ( W) n Q )L/2] (3(oj, n ) , (28) 

and with fio being the center frequency of the SFG spectrum. Assigning £ = n — Oo we then rewrite the photon flux 
amplitude operator (Eq. (|23ll ) as: 

f} SFG {t,L) = =— / d£e- lit g SFG {£) / dw / SPG (w, H ) ai(w)a 2 (n + f - w) . (29) 



2tt 



By applying Eq. (|22p , using ?) SFG as defined in Eq. (|29p . we may evaluate the SFG intensity induced by any initial 
state of the light (note again that this expression is valid only as long as the up-conversion does not deplete the 
incoming fields). 

B. TPA 

Using second-order perturbation theory, very similar expressions can be derived for TPA. The interaction Hamil- 
tonian of an atom and one spatial mode of the electromagnetic field takes the form of: 



j k 

where fikj — (M/J-lj) are the dipole moment matrix elements, Ukj = — <-^j, lk,lj are the level life times, and the 
summation is performed over all the combinations of the unperturbed atomic levels \j), \k). In order to evaluate the 
TPA amplitude, we may use the second-order approximation for the time-evolution operator that corresponds to this 
interaction Hamiltonian, taking only the terms that contribute to a transition from the initial (ground) level \g) to 
the final level |/). Assuming the atom is initially in the ground state \g), the probability Pj FA for a light state \ip) to 
induce TPA can be represented as: 

p/ PA (t) H^ICaW WW (3i) 

with f) TPA (t) defined in a very similar way to Eq. (|29[) : 

oo . rm p i(uj fg -u-u')t 



teW = --, ^y^VfnVng / dw dbj' — ; X (w) & 2 (u/) 

A-Kcebh^ J Q J Q {(jf g -u-u'-i'yf){uj n g-u-i'y n ) 



e -«^ft f e -W [<»> y/u{u, tg + i-- u>) . , ^ A / 

2^M/nPns / " — / dw- — — - — - ai{cj) a 2 {ujf g + t, - w) , (32) 



inceSh ^ J C + *7/ Jo ( w «3 —u — vy n ) 

where the subscripts g, n, f denote the ground, intermediate and final levels, respectively, and £ = lo + lu' — tOf g - Once 
again, the subscripts 1,2 denote the two spatial modes of the fields that induce the interaction. For convenience, 
we assume here that the atom is located at z\ = z 2 — along these modes. For non-resonant TPA, i.e. when the 
spectra of the inducing fields do not overlap with any resonant intermediate levels, the operator defined in Eq. (|32p 
can approximated to: 

AirceSH (ujng - {to}) J £ + ijf J V 
where (uj) is the center frequency of the field in mode 1. 

The similarity between Eq. (f3"3"|) and Eq. (|2"9"|) becomes evident if we denote: 

?)tpa(0 = constant x e'^f"* / d£ e _l?t g T p A (£) / dw /tpa(w, W/ 9 ) ai(w) a 2 (w/ 9 + £ - u>) , (34) 

with: 



o 
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and for non-resonant TPA 



f TFA (u>,UJf g ) — ^ Hf n ^ ng 



(35) 



5tpa(C) 



1 



f TPA (u!,U! fg ) = ^Juj(uj fg -Uj) 



(36) 



where in both expressions we approximated y/u>(u)f g + £ — u>) ps Jw^Wfj — w), since the atomic level linewidth (which 
defines the range over which £ is integrated) is negligible compared to the optical frequencies. 



C. A generic two-photon operator 



In the next section we shall take advantage of the similarity between the expression for SFG and non-resonant TPA, 
and perform all the derivations using the following generic form for the probability amplitude of the final state of the 
nonlinear interaction: 



f)(t,n,n) = e-* > ° t / d£e-^g(0 / du f(u, fi ) a^e^ a 2 (O + £ - u)^ a °+^-^ n . (37) 



By introducing the exponents e _lwl > 2T1 ' 2 we take into account the possibility that the fields di^WjO) have propa- 
gated freely separately, accumulating temporal delays of 71.2, respectively, before inducing the interaction. The final 
probability for the interaction, and hence the intensity I(t, 71,72) of the measured signal is then represented by: 

i(t, n , r 2 ) cx (i/>\ tf(t, n, 7a) v(t, n,T 2 )\i>). (38) 

The effect of any inhomogeneous broadening mechanism of the final level may be taken into account by evaluating 
the intensity Ii^a) of each homogeneously broadened subset with center frequency f^o, and defining: 

jtotai K J p(n )j(fj ) drio 5 (39) 

with P(£lo) being the probability distribution of the center frequency f2o- 

By assigning the appropriate expressions for g(£), f(u>, f2o), Eqs. (|3T|) - ([55)l may represent SFG (using the definitions 
in Eq. (|27|) ), TPA (using Eq. (f35|) ). or non- resonant TPA (using Eq. (f36|) ). as well as other two-photon interactions. 
For example, to evaluate the rate of coincidences of photons at some optical bandwidth A around Oq/2, we may 
assign: 

f(u),Sl ) ps 1 when \tu - Q /2| < A/2 
ps otherwise, 
g(£) PS 1 when |£| < A 

ps otherwise, (40) 

where we assumed that the bandwidth A is smaller than the optical frequency Qq/2. Under these conditions, 
r)(t,Ti,T2) oc Ei (t — TifE^it — 72), and so the overall intensity is simply proportional to the second-order cor- 
relation function between the fields : J(t, t%, 72) cx g^(7i — t%), which is typically taken to represent coincidence events. 

Finally, f(u>, Ho) may also represent any spectral filters Q(u>) that are applied to the inducing light by denoting: 

/(w,fi ) = 0(o;)e(no - «) , (41) 
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where we assume that the spectral filtering is constant within spectral slices that are narrower that the final state 
bandwidth, and so we may neglect the dependence on £. If this is not the case, than we should keep the dependence 
on £: f(ui, Oqj £) = 6(^)6(^o + £ — w ) ■ ln the case of SFG, a spectral filter can be applied to the up-converted light 
as well, in which case its amplitude transmission function should multiply <?(£)■ 

IV. A GENERIC EXPRESSION FOR TWO-PHOTON INTERACTIONS WITH BROADBAND 
DOWN-CONVERTED LIGHT: THE COHERENT AND INCOHERENT SIGNALS 

For the following derivations we shall assume that the signal, idler and pump fields have each a single spatial mode 
and polarization, that the down-conversion process occurred along a distance L, and that the signal and idler may 
then travelled along different paths, resulting in delays of r s ,Ti respectively, before inducing the nonlinear interaction. 
Note that in the following we consider only two-photon interactions that result from cross-mixing of the signal and the 
idler fields and not from self-mixing of the signal with itself or the idler with itself (the possible contribution of such 
self-mixing terms will be considered briefly later on). Assuming that good phase- matching conditions are achieved 
for down-conversion at some bandwidth A DC , we use the expressions obtained for d s ^(oj, L) (Eq. Q17[l) in di^u), 0) of 
Eq. (|37|) . respectively: 

f){t,T s ,Ti) = e- iUot J d£, g(0 e~ lit J du f(u, n^e^' e^+t-^n e -i&k{u>)L 

x I yT+ n s {uj) d s (uj) J— p- \/ n s (u) / duj P A p (u p ) a\(u p - oj) I 

V •'o / 



^n t (n + £-uj) / dw' P k p {uj' p ) al(w' p + oj - fi - £) ■ (42) 



Operating on the initial vacuum state and using a(u))a'(ui') = S(u> — u>') + cv (us')a(u)): 

fj(t,T s ,Ti)\vac) = e~ mot J d£ 5 (0 e" l «* J du f(u, n ) e *-r. e i(n -rt-«)r 4e -iA*(a»)i 

/ _j r°° r°° 

X ( 2^J~ V n s(^) n i(^o + £ ~ w ) J dujp k p {uj p ) J duj'p A p (uj' p ) 

x d\{uj' p + uj — il Q — £)al(w p — uj)\vac) 

+ n s (u))) rii(flo + € - A p (fl + £) \vac) ] . (43) 



The intensity I(t,T s ,Ti) can therefore be separated to two components, which we will denote as the 'coherent' and 
the 'incoherent' signals: 

I(t,T s ,Ti) CX (vac\f]' ! (t,T s ,T l )fj(t,T s ,T i )\vac) = I C (t,T s ,Ti) +P C {t,T Sl T l ) , (44) 

with: 



2irl p 



d£ 9(0 A p (O + e-^-r*) I du f(u, n ) e -* Ak ^ L 

2 







In order obtain P c (t), we shall change variables in Eq. (|43|) : u = lj p — u>, v = u>' + u> — £, leading to: 



(45) 
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e in ( Tt -t) j-oo 

f°(t,T s ,Ti)\vac) = — — - — / (III 



x f(u p - u, n t) )e- lAk ^"-^ L ^n s {Lj p - u) ni (n + £ - uj p + u) 
x A p (u! p ) A p (v + u + £ — uip) a\(v — tlo) K(u) \vac) , 



(46) 



and so: 



I ic (t,r s , n ) 



d£, g{0 e^ (Tl ^ I dujp ^ n s {uj p - u)rii{Q, + £, - uj p + u) 



frhYh du L dv 

x j(up - u, n ) e - lAk{uJ »- u)L A p (wp) A p (v + u + £-ujp) e 1 ^"^-^ 



(47) 



Equations 1451 and 1471 can be calculated numerically by assigning the appropriate /(w, J7o), <7(£) , the spectral ampli- 
tude of the pump A p (w p ) and the power spectrum of the down-converted light n s ^(co). However, these expressions can 
be further simplified analytically by making a few more reasonable assumptions, according to the specific nonlinear 
interaction that is being evaluated. Our only assumption so far was that the pump bandwidth S p is significantly 
smaller than the down-converted spectrum A DC . This assumption enabled us to neglect the dependence of Afc on 
spectrum of the pump, and replace w p , with its center value {lo p ). Similarly, we will assume that n(w) and /(w) are 
also roughly constant within spectral slices that are narrower than the pump bandwidth S p . In order to simplify Eqs. 

and (|4"7]) , we will assume from here on that the bandwidth 7 of the final-state g(£) is also significantly smaller 
than the down-converted bandwidth: 



A DC > 5 P , 7 



(48) 



In TPA 7 represents the bandwidth of the final level 7/, and in SFG it represents the possible (phase-matched) 
bandwidth for the up-conversion process. In accordance with this assumption we will neglect the dependence of n Si i 
on £, which leads to: 



I C (t,T s ,Ti) 



2nl n 



1 



dv 



du 



f({u! p ) - u, fi ) ^n s ((ujp) - u)ni(Cl - (u p ) + u) 



d£ ff(0 e~ l « (t - ri) / dtop A p (lu p ) A p {v + u + £ - lu p ) e 1 ^ 



P c (t,T s , n ) 



(49) 

To further clarify these expressions, let us further restrict ourselves only to SFG and non-resonant TPA, for both of 
which we can approximate the amplitude of spectral function /(w, f2o) to be some average value /o over some spectral 
bandwidth, A, when A is taken to denote only the part of this spectrum that overlaps with the down-converted 
spectrum A DC . We will not neglect, however, the phase of f(ui, f2o); specifically - let us assume that the signal and 
the idler have spectral phases of 9 s (u> s ), 6i(u>i), respectively (for example, due to spectral filters or a pulse-shaper). 
Thus we take f(u } Q ) to be: 



/(w.fio) 



favg X 6 



i[8 3 (u)+6i(ti -u)] 



vhen u> lies within the bandwidth A 







otherwise , 



(50) 



with A < A DC . Let us also focus our attention from here on to two-photon interactions with a final-state frequency 
that is close to the pump frequency (as will soon become evident, outside this regime the coherent signal dies out, 
leaving only the incoherent signal): 



|fi - H)| < A 



(51) 
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We now define an "effective pulse" P e {t) which is the power of a pulse with a power-spectrum that is equal to : 
n e (ui) = (n s (ui) + l)ni(n - ui), and spectral phase that is equal to e *[<?a (f>o-w)] : 



Pe(t) = p- 



F'U V(l + n.(«)) m(n - u) e ilo*(")+Wo-*)) e -i±k(. u )L 



(52) 



where the factor of 2ir originates from the symmetric definition of the Fourier transform as : 

x(t) — F~ {x(a;)} = — 1= J X(oj) e~ luJt duj, and with P being the total power (times 2n) of such a pulse with a constant 
spectral phase (a 'transform- limited' pulse): 



Po 



yj (n a (w) + l)n,(f2 - w) <2w 



A 2 (n 2 + n) , 



(53) 



where n = (n s (uj)) (rij (f2o ~~ Ci; )) is 27r times the average of the photon flux spectral density of the down-converted 
light over the bandwidth A (see Eqs. (fT"4]) - (fT8|) ). Due to the normalization by Pq we obtain: 



< P e (t) < 1 , 



(54) 



where P e (0) = 1 is achieved for the un-shaped (transform- limited) pulse. Finally we also take not of the fact that the 
signals depend only on t — T{ and Tj — r s : I(t, r s , Tj) = /(i — Tj, — r s ). 

Using the definitions of Eqs. (f53|) . (f52|) . I c (t,T s ,Ti), P c (t 7 T s ,Ti) take the following form: 



7 C (i - Ti.Tj - T s ) PS I/, 



(I eg | 



|2 A 2 (n 2 + n) 



J ~y 



i((t-Ti) 



Pe(n - T s ) 



1 2 n 



I \t T ii T i T s) ~ /ai) g I , ., , ■> 



oo 

dtt / dv 



Jo 



(55) 



V. THE RESULTS: THE TEMPORAL AND SPECTRAL BEHAVIORS OF THE COHERENT VS. 

INCOHERENT SIGNALS 



A. General 



Equation [55] already reveals most of the unique features of two-photon interactions induced by broadband 
down-converted light that were mention at the introduction. These features are presented in Figs. WSM In order to 
present a quantitative picture of the behavior and the relative magnitudes of the coherent and incoherent signals, 
we assume in the following realistic physical parameters, that are similar to the experimental parameters in [33. l40l|. 
Specifically, we assume broadband, degenerate but non-collinear down-converted light at a bandwidth of 80nm 
around 1033nm, and consider two-photon interactions around f2 = 516.5 nm. Specifically, for TPA we assume a 
final- level bandwidth of 7 TPA rs 5MHz, and for SFG we assume a phase-matched bandwidth for up-conversion of 
7uc = 0.3nm. We assume the down-converted light was generated by a pump that is a ~ 3 ns pulse with a bandwidth 
of 5 p = O.Olnm around 516. 5nm. Such parameters are typical for Q-switched laser systems. 

The first unique feature is the coherent signal's non-classical linear intensity dependence: P cx n 2 + n. This 
behavior manifests the fact that the signal and idler modes share a single wavefunction. Figure [T] depicts this 
behavior of the coherent signal, compared to the quadratic intensity dependence of the incoherent signal. The relative 
magnitudes of the incoherent signals for SFG, TPA and coincidence events are calculated using Eqs. (f6"5| . ([57]) and 
(|96l) . respectively (derived in the following subsections). 

Note that while the dependence of the coherent signal on the flux of the down-converted photons may be linear, 
the response of the two-photon signal (TPA, SFG or coincidences) to attenuation of the down-converted light by 
linear losses (namely absorption or scattering, for example by optical filters or beam splitters) is always quadratic, as 
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FIG. 1: Power-dependence of the coherent and incoherent parts of two-photon interactions with down-converted light, (a) The 
relative magnitudes of the coherent signal I c (black line) vs. the incoherent signals for TPA, SFG and coincidence events (I^ A , 
Isfq and Rc lc , respectively, gray lines), assuming typical physical parameters described in the text, represented on a log-log 
scale. The graph shows the n 2 + n dependence of I c , with n being the average spectral photon density, i.e. the total photon 
flux in each of the signal and idler modes, divided by their bandwidth A. All the incoherent signals demonstrate a quadratic 
dependence on n. The dashed line represents a completely linear function, for comparison, (b) The dependence of I c on n 
represented on a linear scale, demonstrating the nearly linear dependence at small values of n. The dash-dot line depicts the 
quadratic response of I c assuming that down-converted light with n = 0.2 is being attenuated by linear losses. 



is evident from the presence of the term |/ a u S | in the expressions for both the coherent and incoherent signals. This 
behavior is depicted by the dash-dot line in Fig. [TJ which assumes that down-converted light with average spectral 
photon density of n = 0.2 is being attenuated by optical filters. These results are in excellent agreement with the 
experimental results of SFG with entangled photons presented in [24[ . 

The second unique feature that appears in Eq. {55)) is the pulse-like response of the coherent interaction to a relative 
delay between the signal and the idler, as represented by the term P e (r; L — t s ), which is the (normalized) response 
of mixing two ultrashort pulses with practically the same spectra as the signal and idler. As such, P e (ji — r s )) is 
sensitive to dispersion, including the dispersion that was accumulated in the down-conversion process itself, denoted 
by the term e - lA H") L m the definition of P e (ji — t s ) (Eq. (f52"| ). Figure[^c) shows a zoomed-in picture of this sharp 
temporal dependence of the coherent signal on relative delay between the signal and the idler beams, assuming that 
dispersion is either negligent or is compensated by spectral phase filters, leading to a sharp response which is exactly 
as if the interaction was induced by a pair of 35/s (transform-limited) pulses. 

The coherent summation over the spectrum which leads to this pulse-like behavior also implies that the coherent 
signal can be shaped by spectral-phase manipulations, exactly like with coherent ultrashort pulses. Note that the 
shape of P e {ri — t s ) is determined by the sum of the phases applied to antisymmetric spectral components of the 
signal and the idler: 9 s (ui) + 6*;(f2o — <jj)- This implies that if the same phase filter is applied to both the signal 
and the idler beams (or the same dispersive medium), only spectral phase functions that are symmetric about 
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FIG. 2: The temporal behavior of the coherent and incoherent parts of two-photon interactions with down-converted light, as a 
function of a relative delay between the signal and the idler fields, (a) The relative magnitudes of the coherent signal I c (black 
line) vs. the incoherent signals for TPA, SFG and coincidence events (Jtpa, Isfg an d Rc lc , respectively, gray lines), assuming 
typical physical parameters (described in the text), represented on a logarithmic scale. This graph depicts the instantaneous 
peak of the signals at t = n for one typical down-converted pulse, and shows the contrast between the sharp temporal behavior 
of I c , and the behavior of the incoherent signals, which is on the same ns-timescale as the pump pulse. Note that Itpa presents 
a smooth and symmetric graph, reflecting fact that the long-lived (~ 30 ns) final atomic state integrates over the intensity 
variations of the pump pulse. Similarly, Rc c , which takes into account integration over the gating time Tg ~ 1 ns is more 
smooth than Ism, which represents an instantaneous parametric process. For the case assumed here of a non transform-limited 
3ns pump pulse, the ensemble average of all the incoherent signals becomes proportional to the normalized intensity-correlation 
function of the pump g^ 2 ' (r< — T a ), which exhibits a ~ x2 'bunching peak' at delays which are shorter than the pump's coherence 
time 2n/S p . This is shown in (b) in black, together with the instantaneous incoherent TPA signal I^ A (gray), which can be 
approximated to be proportional to g^ 2 ' even without averaging, (c) A zoomed-in presentation of the ultrashort-pulse like 
behavior of the coherent signal: I c (ji — r s ) oc P e (ri — r s ), assuming dispersion is compensated and no additional spectral 
filtering, resulting in a behavior that is identical to that of mixing two 35 fs transform-limited pulses, (e) I c at zero signal-idler 
delay as a function of the magnitude of a spectral phase filter 6 s {uj) that is applied to the signal, for example by a pulse-shaper. 
s (ui) is drawn in (d) (black line) at a magnitude of tt, together with the down-converted power spectrum (gray line), (f) The 
shaped temporal behavior of the coherent signal with 8 a (uj) as depicted in (d) applied to the signal field. 



fio/2 ~ (w P ) /2 would affect P e (ri — r s ). Figures [Hc)-(e) demonstrate how applying a spectral phase filter to the signal 
(or idler) spectrum leads to the same result as with coherent ultrashort pulses. This ultrashort-pulse-like behavior 
(including the ability to tailor it by a pulse-shaper) was demonstrated experimentally with high-power SFG in [3^ |. 
high power TPA in [4(j, and with broadband entangled photons in [4l[, with excellent agreement with our calculations. 

Interestingly, P e (ri — t s ), does not depend on the specific type of the two-photon interaction, i.e. the coherent 
signal of TPA, SFG or coincidence events will always exhibit this ultrashort-pulse like behavior. The contrast 
between the temporal behavior of the coherent signal and that of the incoherent ones is shown in Fig. [21a). Note 
that Fig. [DJa) presents the instantaneous peaks of the coherent and incoherent signals at t = Ti for one, single-shot 
arbitrary example. As is evident, the incoherent signals (calculated for TPA, SFG and coincidence events by Eqs. 
(ItTTI) . (|T5j) and (1§5"|) , respectively, assuming n ^> 1) always demonstrate a temporal dependence on r s — n that is on 
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the same ns-timescale as the pump pulse. 

It is important to note that the duration of such Q-switched pulses is much longer than their coherence time (which 
in this case is 2n/S p sa 89 ps), which was their duration if they were transform- limited. Such pump pulses, for which: 



t p > 2ir/S p , (56) 

can be considered a 'quasi-continuous' light, since they can be viewed as short bursts of continuous light, especially 
when time-scales that are shorter than r p are considered, during which the average intensity stays roughly constant. 
Thus, such quasi-continuous pump pulses yield approximately the same results for I c , I lc as a continuous pump, 
especially when the ensemble average of many such pulses is considered. In particular, as will be shown in the next 
subsections, once averaged the incoherent signals all becomes proportional to the normalized second-order correlation 
function of the pump gp (t, — t s ). This behavior is depicted in Fig. [DJb), which depicts the calculated gp 2 '(Ti — t s ) 
together with a zoomed- in presentation of the instantaneous incoherent TPA signal /^pa- As i s evident, even without 

(2) 

averaging I*p A follows very closely g p (r, — t s ). This is explained by the fact that the incoherent excitation of the 
long-lived (~ 30 ns) final atomic state actually averages the intensity fluctuations of the pump (see Eq. (|79| V When 
the ensemble average of many such quasi-continuous pulses is taken, the incoherent TPA signal, as well as and 

( 2 ) 

Rc lc become practically identical to g p '(ji — t s ), demonstrating the expected x2 'bunching-pcak' at delays which 
are shorter than the coherence length of the pump (r s — Tj < 1/S P ) [4j. 

Another result of the temporal integration which is performed by the incoherent TPA process, is the fact 
that the behavior of I^ A (Ti — r s ) is more smooth and symmetric than that of X^fg — r s)> which represents an 
instantaneous process. Since coincidence measurement also includes a temporal integration over the ~ 1 ns gating- 
time of the detectors, Rd c presents a behavior which is more smooth and symmetric than I^ G , but not as much as I^pa- 

Unlike the incoherent signals, the coherent signal's sharp behavior P e {ji — t s ) depends on the large-scale properties 
of the down-converted spectrum, and therefore its shape is not affected by shot-to-shot noise, only its relative height. 
Thus, we see that there are three timescales in our system. One is the duration of the pump pulse (which can be 
infinity for a continuous pump), the other is the coherence time of the pump which is ~ l/$p (and is equal to the 
duration of the pump pulse, in case it is a transform-limited one), and the shortest time scale is the behavior of 
the coherent signal, which is on the same timescale as the coherence time of the down-converted light: 1/A. In the 
case considered in Fig. ® the 35/s pulse-like behavior of the coherent signal stands in contrast with the temporal 
behavior of the down-converted light itself, which is a 3ns pulse in this case, i.e. 85,000 times longer. The effect is of 
course even more intriguing when continuously-pumped down-conversion is considered. 

The temporally-sharp behavior of the coherent signal also stands in contrast with its equally sharp spectral behavior. 
More specifically, we need to distinguish between two spectral behaviors. One is the excitation-spectrum, i.e. the 
frequencies which are excited by the two-photon interaction. The other is the dependence of the interaction on the 
pump wavelength. In the case of SFG, the excitation spectrum corresponds to the spectrum of the up-converted 
light. In the case of TPA this spectrum corresponds to which atomic levels will be excited. By Fourier-transforming 
the amplitude of the coherent signal in Eq. ([55]) back to the frequency domain O of the generated signal, we see 
immediately that the excitation power-spectrum of the coherent signal is simply the spectral overlap between the 
narrowband pump and the final state, and does not reflect the broad spectra of the signal and the idler fields which 
induce the interaction: 



, |2 A 2 (n 2 +n) 

| Ja-vg | T *e \ 'i 



ff(O-Oo) A p (ft)e 



(57) 



This implies that if the pump bandwidth is narrower than the final state bandwidth, the excitation spectrum 
I c (£l) would follow that of the narrowband pump, as is shown for SFG in Fig. [3j In other words, the coherent signal 
behaves as if the pump itself was inducing the interaction. While the spectral behavior of the incoherent signal is 
harder to deduce out of Eq. (|55[) . in the following we show that it is approximately that of the final state; this is 
shown more easily if we assume the final state is significantly broader than the pump (Eq. (|61[) ). or the other way 
around (Eq. (|74[) ). This spectral behavior of the coherent and incoherent was demonstrated experimentally in (38l. l39j| . 



For the case of TPA, even if the pump is narrower than the final atomic state, this is not reflected in the spectrum 
of the fluorescence from that level, because the temporally random, incoherent emission process of the emission erases 
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FIG. 3: The spectrum of the coherent and incoherent parts of SFG light induced by high-power (n ^> 1) broadband down- 
converted light, with physical parameters as described in the text. The graphs are shown on a logarithmic scale in (a) and a 
linear scale in (b), and demonstrate how the coherent signal (black) replicates the narrow spectrum of the pump (S p = O.Olnm). 
In contrast, the spectrum of the incoherent signal (gray) is as wide as the phase-matching conditions at the crystal allow it to 
be (7u C = 0.3nra). 



the information on the exact frequency that drove the transition (especially in this limit of weak, non-stimulated 
interaction). However, while the excitation spectrum may not be directly accessible, the other kind of spectral 
behavior, i.e. the dependence of TPA on the pump wavelength, can be explored experimentally. As already evident 
from Eq. (|55[) . the magnitude of the coherent signal depends on the total spectral overlap between the pump 
spectrum and the final state spectrum. This is shown in Fig. 2J which considers TPA in atomic Rubidium (Rb) at 
the 55*1/2 — 41)3/2,5/2 transition. The coherent TPA rate indeed behaves as if the pump laser itself was inducing the 
transition (which is of course a forbidden transition for a one-photon process), demonstrating a spectral resolution 
of O.Olnm, almost resolving the IZAGHz hyperfine splitting between the 4D 3 / 2 and the 4D 5 / 2 levels, even though 
the interaction is induced by a light with a total bandwidth that is ~ 2000 times wider. For the sake of simplicity 
we ignored in Fig. [4] the hyperfine splitting of the 55 , i/ 2 ground-level in Rb (3GHz for Rb87, 6.8GHz for Rb85). 
This calculation should be compared with the experimental results of [40] (in which a wider pump bandwidth of 
0.04nm and power-broadening prevented resolving the hyper-fine splitting, nonetheless demonstrating a spectral 
resolution that was ~ 2000 times narrower than the down-converted bandwidth). In contrast, the incoherent signal 
(for SFG, TPA and of course coincidence events) is practically independent of the wavelength of the pump (see Eq. 
(frTTj) for SFG and Eq. (f74")) for TPA). The incoherent signal responds only to the change in the down-converted 
power spectrum that results from the change of the pump wavelength, and so it exhibits only the very wide spectral 
response that is expected in an interaction that is induced by mixing of two 80nm wide incoherent beams. This can 
be viewed as resulting from the fact that the incoherent interaction has no knowledge of what was the wavelength of 
the pump that generated the down converted light, since that "information" lies only in the correlations between the 
spectral phases of the down-converted modes - phases that play no role in the generation of the incoherent signal. 

In the following subsections we present more specific analytic expressions for the temporal behavior of the coherent 
and incoherent signals, performing approximations that fit various two-photon interactions. First we consider the 
case where the spectral width of the final state exceeds that of the pump, as is typically the case with SFG and is 



15 



1 

0.8 

0.4 
0.2 


516.6 516.62 516.64 516.66 516.68 516.7 
Pump center wavelength [nm] 

FIG. 4: The dependence of the coherent part of TPA (I^ PA , black) induced by down-converted light, as a function of the center 
frequency of the pump. The graph shows how I% PA behaves as if the O.Olnm-wide pump itself is inducing the transition (5Si/ 2 
— > 4D 3 /2,4D^/2 in atomic Rb, with the transition wavelengths drawn in gray). The incoherent part of the TPA is practically 
insensitive to the pump wavelength, and since for practical power levels (n > 1) it is smaller than the coherent part by roughly 
8 P /A » 1/2000 (with 8 P , A being the bandwidths of the pump and the down-converted light, respectively), it is too small to 
appear in this graph. 



always the case with coincidence detection. Then we consider the case where the spectral width of the final state 
is smaller than that of the pump, as can occur with TPA. For more accurate results, and for experimental schemes 
that do not comply completely with the assumptions and approximations that follow, Eqs. (HH) or ([53)1 should be used. 



B. Pump bandwidth smaller than the final-state bandwidth (example: SFG) 

So far we have assumed that the down-converted bandwidth is significantly larger than the pump bandwidth and 
the final-state bandwidth (Eq. (J48J) ) . In the following derivation we consider the case where the pump bandwidth is 
also significantly narrower than the final-state bandwidth: 

A > 7 > S p . (58) 

This condition allows us to simplify the expressions for both the coherent and the incoherent signals from Eq. (|55|) 
by using the following approximation: 

/ d£g(Q A P (cD + C) e-^-^ « g((w P ) - w) A p (t - n) e io ^ , (59) 
where A p (t) = F~ {A p (lj)\ is the temporal amplitude of the pump, leading quite immediately to: 

r(t-n,Ti -t b ) 

r c (t-T u n-T s ) 

with I p (t) — \A p (t)\ 2 . To represent more accurately typical experimental conditions, we can take the ensemble 
average of the signals (i.e. averaging over many pulses in the case of a pulsed pump, or over time in the case of a 
continuous pump). In particular, if we consider a continuous pump or a quasi-continuous pump (i.e. the center part 
of a non-transform- limited pump pulse for which t p 3> 2ir/6 p ), we can approximate for t — r Sl t — r, <C t p : 
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(60) 
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(I p (t - n)) « <I P (« - r,)) « <I P (0)) « / p . (62) 
This allows us to further simplify the expressions for the coherent and incoherent signals by using: 



(63) 



(i p (t - n)i p (t - t s )) (i p (t - n )i p (t - t s )) 

fi ~ 7T77 u/t 71 yT = Sp( T » _ r «) ' 64 ) 

where g 2 (ri — r s ) is the normalized second-order correlation function of the pump field. Note that the pump field was 
taken as a classical amplitude throughout this paper, hence g 2 represents in our derivation only classical intensity 
correlations. This leads to: 

(I c (t - n, Tl - T S )) w \g((u p ) - n )\ 2 \f ava \ 2 A 2 {n 2 + n) P e ( Tl - t s ) . (65) 
(P^t-T^n-Ts)) « [ d£\g(£)\ 2 \f avg \\ 2 p (T t -T s )An 2 . (66) 

As is reflected in the term J d£' | g(£') | 2 in Eqs. (|6Tj) , (|66|) , the power of the incoherent signal is indeed the incoherent 
summation of all its spectral components, and its spectrum follows that of the final state. This behavior is depicted 
in Fig. [3] for SFG, in which case we may write: 

\favgf \g{Q - n )\ 2 = L 2 smc 2 [Ak((Lu),n)L/2] f3 2 {(u),Q) . (67) 

With the exception of very long crystals (i.e. more than ~ 1 cm in the case of type-I up conversion, or more than 
a few mm in the case of type-II up conversion), the up-converted bandwidth 7 UC is typically at least of the order of 
0.01 THz, and could even reach tens of THz for very short crystals. This bandwidth is orders of magnitude larger 
than the typical bandwidth of continuous lasers, and even larger than the bandwidth of Q-switched, ns-pulsed lasers. 
Thus, the condition of Eq. ([55)) is typically satisfied for SFG, unless ultrashort pulses (hundreds of fs or less) are used. 

Equations 1611661 displays both the spectral and the temporal behaviors of the incoherent signal more clearly than 
Eq. (f55|) did. As is evident, the incoherent signal is practically insensitive to the pump wavelength, and its response 
to a relative delay between the signal and the idler is very slow, since it depends on the temporal behavior of the 
pump, which is either continuous or a very long pulse, compared to the ultrashort-pulse-like behavior of the coherent 
signal; as demonstrated by the term I p (t — Ti)l p (t — r s ), the incoherent signal depends only on the temporal overlap 
between the intensities of the signal and the idler (which follow the intensity of the pump). This behavior is depicted 
in Fig. [5](a)-(c). In particular, Fig. [2Kb) shows the 'bunching' peak which is exhibited by the incoherent signal for a 
quasi-continuous pump at delays which are shorter than the pump's coherence length. 

Since the observed signal is always the sum of the coherent and incoherent contributions, it interesting to compare 
the magnitude of the coherent signal to that of the incoherent one. The average ratio between the coherent signal 
and the incoherent one is therefore: 



(Pjt-T^n-Ts)) (l p (t-Tj))l p A jn 2 + n) 

(l ic (t-n,Ti-T s )) ~ (Ip(i - n)lp(f - r s )> 7 n 2 

which for a continuous or quasi-continuous pump becomes: 



Pe(n -t s ) , (68) 



(nt-run-T.)) A (n 2 + n) 

(Iic(t-T h n-T s )) 7g 2 (r,-r s ) n 2 ^ ^ " (&9j 

In the absence of spectral phase filters, and assuming that dispersion is corrected, we can assign P e (0) = 1- Thus, 
taking into account that typically 1 < S P ( T i ~ T s) < 2, the expected ratio between the coherent and incoherent 
signals approaches A/7 for high powers (n > 1), and even more for low powers where n + n 2 3> n 2 . Thus, under 
the conditions of broadband down conversion and a narrowband final state assumed throughout our derivation, the 
coherent signal dominates over the incoherent signals (see the relative magnitude of I]£ a in Figs. [1][3]). 



17 



C. Pump bandwidth larger than the final-state bandwidth (example: TPA) 

For the following derivation we assume that the final-state bandwidth 7 is significantly narrower than the pump 
bandwidth, as may often be the case with TPA: 



A > S p > 7 . (70) 

Assuming this, we can simplify the expressions from Eq. (|55[) by assuming that the spectral amplitude of the pump 
remains constant within spectral slices which are narrower than 7 : 



(71) 



Additionally, in order to evaluate the incoherent TPA signal we use Parseval's theorem to obtain the following 
relation: 



dv 



dcop Ap(wp) A p (v - Lu p )e ZUJpT 



/oo 
Ip(t' +T)I p (t') dt' 
"OO 



leading to: 



(72) 



rcf, \ \t I 2 I a /o M 2 A2 (n 2 + n ) 

I C {t-Ti,Ti-T s ) « \f avg \ \Ap(n )\ " 

ZlTln 



dHg(0 e 



-iS(t-n) 



Pe(n - T s ) 



(73) 



,9 An 2 f°° r 
pc(t- Tl ,T t -T s ) ^ \f avg \ 2 — J I p (t> + r)lp(t') dt> Jc%g(0 



(74) 



Equations l73ll74l show that in this case the spectra of both the coherent and the incoherent signals are determined 
by the spectrum of the final state <?(£). 

To clarify the temporal behavior of both signals, let us use 



9(0 e-*^ = G(t - n ) , (75) 
with G(t) — F~ 1 {g({;)} being the slowly varying envelope of the temporal response of the final state. This leads to: 



i c (t - n,n — T s) ~ \favg\ |A p (n ) 



2 A 2 (n 2 + n) 



G(t-Ti)\ P e {n-T s ) 



(76) 



,9 A n 2 r°° ,9 
I lc (t-T tl T 2 ~T s ) » \f avg [ — / I p (t' + r)Ip(t')dt' |G(t-r<)r 

p J — OO 



(77) 



Considering TPA as a probable example for the case where the final state is considerably narrower than the pump, 
we can substitute = . Using F~ 1 {g({;)} = v(t — Tj) e iM* - Ti > (where u(t) is a step- function), we obtain for 
TPA: 



r(t-Ti,Ti ~T S ) = f% 



|Ap(fio) 



u(t - n) e - 2 7/(*-0 A 2 (n 2 -(- n) P (t, - r s ) 



P c (t - r l; r, - r s ) = Ktpa * l} ' u(t - n) e - 2 7/(*-^) A , 



I 2 

p 



(78) 
(79) 
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with 



Ktpa — 



32tt 3 



ceSh 



(uHVfg ~ (")) 
(w„ 9 - (uj)) 2 



(80) 



In order to have a better estimation of the magnitude of the coherent TPA signal as compared to the incoherent one, 

I 1 2 

let us try to estimate the average magnitude of A p (17o) , assuming of course that f2o is within the pump spectrum 
(i.e. |f2o — (^p)l < Sp/2). For a pulsed pump, we can use Parseval's theorem again to relate between the average 
spectral power of the pump ^|A p | 2 ) to the average temporal power I p : 



5 P (|A P | 2 > = T P I P . (81) 

However, to get a similar relation for a continuous pump we need to identify the longest timescale in our system, 
which defines the smallest frequency increment, i.e. the quantization unit of our frequency domain. In the case 
considered in this subsection, the smallest frequency scale is that of the final state, 7/ (Eq. ([70])). hence the longest 
relevant timescale is the final state lifetime 77. Accordingly we can approximate for a continuous pump: 



MN 2 ) « T flp, (82) 

which is identical to Eq. (|81[) . only with 77 replacing t p . In other words, although the pump is continuous, we can 
treat it as if it was composed of serious of pulses, each of them 77 long. Since the atomic state, which is the slowest 
component in our system, has a "memory" only 77 long, its response is not affected by interactions that occurred 
more than 77 seconds ago. All the other components of our system have shorter coherence time; for example, since 
the pulse bandwidth is significantly wider than 7/ = 2w/Tf, it is not affected by such 'chopping' of continuous light 
to pulses Tf long. Using Eq. (|82[) . and assuming the final state is included within the pump spectrum, we get: 

I c (t - Tl , Tl - T .) « k tpa ^ u(t - n) e" 2 ^*-^ A 2 (n 2 + n) P e (n - r s ) , (83) 

dp 

and the ratio between the coherent and the incoherent signals is therefore: 



I C (t ~ Tj,Tj - t s ) ^ Ip A 

Pc(t-T U T % -T s ) ~ ^Ipit' + T S - n )lp{t') dt' 5p 



n 



Pe(n-T s ). (84) 



Similarly to the case with SFG, we see that the coherent signal is stronger than the incoherent one (in the absence 
of dispersion or delay between the signal and idler), this time roughly by the ratio between the down-converted 
bandwidth and the pump bandwidth A/S p . 

To clarify the dependence of the incoherent signal on a relative delay between the signal and the idler, let us consider 
the case of a continuous, stationary pump, for which (l p (£ + r) ) = (lp(t)) = I p , and so the intensity correlations can 
be represented by the normalized second-order correlation function of the pump: 



f ip(t' + T s -n)ip(t>)dt> (ip(t + T s -n)ip(t)) 

~ Tl ^ ^ (i p(t + r s -n)>(i P (0> = Tf 4 ( " " Ts) ' (85) 

where 77 is the lifetime of the atomic |/) state, which is the physical time interval over which the intensity correlations 
are actually integrated. Equation [85] is valid as long as the lifetime of the final state is much longer than the coherence 
time of the pump, which is indeed the case considered here since we assumed S p ^> 7. Therefore we may write: 

r c (t - Ti ,n - t s ) « Ktpa Tf gp 2 Hn - t s ) u(t - n ) e - 2 ^*-^ A n 2 , (86) 

and the ratio between the coherent and incoherent signals then becomes: 
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If the final state is inhomogeneously broadened with an inhomogeneous bandwidth 7 IH , we may use Eq. (l39l) to 
obtain (assuming this time that 7lH > 5 p ): 

I^(t - n,n - t s ) a k tpa Tf u(i - n) e^C*-*) A 2 (n 2 + n) P e (T< - t s ) 

I£(t - n, n - r.) « Ktpa r y gf fa - r s ) u(* - tj) e^*"^ 7lH A n 2 , (88) 
and the ratio between the coherent and the incoherent signals is then the same as in Eq. (|87j) . only with 7lH replacing 6 P . 

Similar results can be obtained for a quasi-continuous pump, i.e. if we consider the ensemble average of non 
transform-limited pulses, for which we may write (for small delays, Tj — t s <C t p ): 

( ^ M ' ,+ y )M ' v > ^>(,-,). m 

This leads to the same results as with a continuous pump, with the only difference being that t p is replacing 77: 



(r(t- n ,n-T s )) « KTPA ^u(i-T 4 )e- 2 ^^)A 2 (n 2 +n)P e (r l -r s ) (90) 

Op 

{P c (t-n,n-T s )) « KTPATpgftn-r^u^-r^-^^An 2 , (91) 

and the ratio between the average coherent and the average incoherent signals therefore remains the same as in 1871 
Accordingly, in the case of inhomogeneous broadening the same results hold, with t p replacing Tf there as well. 

D. Coincidence events 

It is quite intriguing to compare the results obtained for SFG and TPA with down-converted light to the expected 
rate Rc of coincidence events, i.e. the simultaneous arrival of signal and idler photons. Typically, the coincidence 
rate is evaluated as proportional to the second-order correlation function g^(i"). If the temporal response of the 
coincidence detectors (and the corresponding electronics) is slower than the coherence time of the photons (~ 1/A), 
as is typically the case, then it is taken into account by integrating over the gating time Tg : 

Rcoc / g (2) (r) dr . (92) 

JTg 

In order to obtain an approximated expression for the coincidence rate, we will use the spectral functions 
<7(£), f(oj, Slo) defined in Eq. (|40|) . Essentially, this means that we treat coincidence detection as if it was an SFG 
process with a very large up-converted bandwidth 7 UC : 

7uc = 2 A . (93) 

Intuitively speaking, such SFG process may be considered as equivalent to coincidence detection since any pair of 
photons that arrives at the crystal simultaneously (i.e. with a temporal separation that is smaller than their coherence 
time ~ 1/A) has an equal probability to be up-converted, regardless of the frequency of the resulting up-converted 
photon. Although we have previously assumed that 7 UC <C A (Eq. (|48p ). this assumption was made only to allow the 
neglect of the spectral variations of n Sj i(u>). Therefore, if we limit our discussion to down-converted spectrum that its 
average is approximately smooth, we may use the expressions obtained for SFG to describe g^ 2 ^(r s — t»), simply by 
replacing 7 U0 by 2 A: 

^(r.-rO oc A 2 ^ + n)P e { Tl -T s ) 

where the first term represents the coherent contribution, and the second represents the incoherent one. However, 
g( 2 '(T s — Ti) represents the actual coincidence detection rate only for infinitely fast detectors; for broadband radiation 
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the temporal resolution of the coincidence detectors is typically orders of magnitude longer than the coherence time 
of the photons. Therefore, applying Eq. (J^U) and taking Tg ^> 1/A, \t s — Tj|, we obtain: 

Rc c K Mt-n) A(n 2 + n)Pe(T ._ Ts) 

Rc iC OC ^ft-^M*-^) 2rgA 2 n 2 _ (Q5) 
V 

As with TPA or SFG if a temporal averaging is performed for the case of a continuous pump, or an ensemble average 
is performed with quasi-continuous, long pump pulses, the ratio (l p (t — Ti)/I p \ in the coherent term approaches 1, 
and the ratio (l p (t — Ti)l p (t — t s ) / 1%) in the incoherent term approaches gp(r,; — t s ) (for small delays). Thus the 
ratio between the average coherent and the incoherent contributions becomes: 

Rc c 1 1 (n 2 +n)P e (n-T s ) 

Rc lc ~ gp( r i - t s ) 2Tg A n 2 1 j 

Since at any power level down-converted light is essentially composed of simultaneously-created photon pairs, it is 
natural to assume that it will exhibit high degree of bunching, in the sense that there will always be a significantly 
higher rate of simultaneous arrivals of photons from the signal and the idler beams (at zero delay), as compared to 
Poissonian or even thermal distributions. However (and counter-intuitively), as is evident from Eqs. (|95|) - ([96l) , since 
Tg ^> 1/A, at high power levels (n > 1) the coincidence rate is dominated by the incoherent term, which exhibits 
similar bunching properties as those of the pump; thus, unlike the frequency-selective processes of SFG and TPA, 
the coherent contribution to the coincidence rate is dominant only at the very low power levels of n <C 1, where 
down-converted light can be described as a stream of entangled photon pairs. 



VI. SUMMARY AND CONCLUSIONS 



In this paper we derived expressions for two-photon interactions induced by broadband down-converted light 
that was pumped by a narrowband laser. In section [II] we solved the equations of motion for the annihilation and 
creation operators of broadband down-converted light generated by an arbitrary narrowband pump. In section Mil 
we formulated operators that represent the photon flux or the probability amplitude of weak two-photon interactions 
(i.e. assuming low efficiency of the interaction, so that the inducing fields are not depleted) induced by arbitrary 
broadband light. In sect ion [TVl we combined the results of the previous sections to obtain expressions for the intensities 
of two-photon interactions, namely SFG, TPA and coincidence events, induced by broadband down-converted light, 
and in section [V] we explored their temporal and spectral behaviors under various conditions. 

Our calculations show that the intensity of two-photon interactions induced by broadband down-converted light 
can be represented as the sum of two terms, one (P) that exhibits a coherent behavior, and a second one (P c ) that 
exhibits an incoherent behavior: 

jtotal _ jc _|_ jic _ ^gy-j 

The two terms vary dramatically both in their spectral properties, as well as in their temporal properties. We 
considered the case where the signal and the idler may propagate freely along different optical paths from the down- 
converting crystal, accumulating independent temporal delays r s , n before inducing the two-photon interaction. The 
coherent signal then responds to a relative delay between the signal and the idler in an ultrashort-pulse like behavior: 

P(t, T „Ti) « i p (t - n ) P e (n - t s ) , (98) 

where I p (t) is the power of the pump (in units of photon flux), and P e (ji — t s ) is the temporal response one 
would have got if the two-photon interaction was induced by mixing two ultrashort, transform- limited pulses with 
the same power spectra as the signal and the idler, although the signal and the idler are each incoherent and 
may even be continuous (see Fig. (He)). Accordingly, P e {n — t s ) is sensitive to dispersion just as a coherent 
ultrashort pulse (including the dispersion of the down-conversion process itself), and can even be shaped by 
conventional pulse-shaping techniques (see Figs. [U[d)-(f)). Note that at low powers this corresponds to shaping 
of the second-order correlation function g( 2 ) of the down-converted entangled photon pairs. It is also interesting 
to note that P e (ji — t s ) responds to the antisymmetric sum of the phases applied to the signal and the idler: 
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6 s (oj) + 9i((oj p ) — w), with (u p ) being the center frequency of the pump. Thus, if the same spectral filter is applied 
to both the signal and the idler beams, only phase functions that are symmetric about (u) p )/2 affect P e . Simi- 
larly, if the signal and the idler travel through the same medium, only odd orders of dispersion will have an effect on P e . 

In contrast, the incoherent signal depends only on the temporal overlap between the intensities of the down- 
converted signal and idler beams, and so reacts to a delay between the signal and the idler on the same time-scale as 
the long pulses or even continuous behavior of the pump: 



P C {t,T s ,Ti) CX I p {t - Ti)l p (t - T S ) (99) 

if the pump is narrower than the final state, or 

n OO 

P C (t, T slTi ) CX / ijf + T S - n)i p {t') dt' , (100) 



— oo 



if the final state is narrower than the pump. In both cases, if we consider the temporal average in the case of a 
continuous pump, or the ensemble average in the case of a quasi-continuous pump (i.e. non-transform limited pulses, 
for which r p 3> 2tt/5 p , with r p being the duration of the pulses, and S p their bandwidth), then the average incoherent 
signal is proportional to the normalized second-order correlation function of the pump (for Tj — t s -C t p ): 

{P C (t,T Sl T t ))^g 2 p (n-T s ) . (101) 

Thus, as depicted in Fig. [2j there are three temporal timescales in our system. The longest one is the duration 
of the pump pulse (which can be infinity for a continuous pump). This timescale dictates the temporal behavior of 
the incoherent signal as a function of the signal-idler delay. The next is the coherence time of the pump which is 
~ 1/S P (and is equal to the duration of the pump pulse, in case it is a transform- limited one). For signal- idler delays 
which arc shorter than this coherence time, the average of the incoherent signal is higher since the intensities of the 
signal and the idler become correlated, as they both reflect the intensity fluctuations of the pump. The shortest time 
scale is the behavior of the coherent signal, which is on the same timescale as the coherence time of the broadband 
down-converted light: 1/A. 

As for the spectral behavior, the coherent signal behaves as though the interaction is actually being induced by the 
pump itself, and not by the down-converted light. Thus, the coherent signal is induced only if the pump spectrum 
overlaps with the final state: 



i c (n + « 



f d£ 5 (0A p (f> + 



(102) 



where g(£) represents the the spectrum of the final atomic level in TPA, or the phase-matching function for up- 
conversion in the case of SFG, and with Qq being the center frequency of the final atomic level or of the phase-matched 
spectrum in case of SFG. The consequences of this spectral behavior is that by scanning the pump wavelength we can 
perform two-photon spectroscopy with the spectral resolution of the narrowband pump, even though the interaction 
is induced by light that is orders of magnitude wider than the pump, and not by the pump itself (see Fig. 2]). In the 
case of SFG this means that light is being up converted only at those wavelengths: 

^ FG (f!o + 0«|s(OA P (fto+0| 2 , (103) 

so that even if the phase matching conditions allow broadband up-conversion, Jg FG replicates the narrow spectrum of 
the pump (see Fig. [3]). 

The incoherent signal, on the other hand, is insensitive to the exact wavelength of the pump that generated the 
down-converted light. Since the information on the original wavelength of the pump is imprinted in the phase 
correlations between the down-converted modes, it affects only the coherent signal P. Accordingly, the incoherent 
signal is induced at all the possible frequency band of the final state of the interaction (see Fig. [3j: 

P c (n + 0^\g(0\ 2 ■ (104) 

The coherent and incoherent signals also exhibit different dependencies on n, the average photon-flux spectral 
density, and on the bandwidth A of the down-converted light. While the incoherent signal depends quadratically on 
n, the coherent signal includes an additional, non-classical term that depends linearly on n: 

P cx (n 2 +n) (105) 
P c cx n 2 . (106) 
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This behavior is presented in Fig. [TJ 

Additionally, since the coherent signal results from coherent summation over the entire (correlated) spectra of the 
signal and the idler, it depends quadratically on A, while the incoherent signal depends only linearly on A 

Thus, excluding the case of pump pulses that are transform-limited, the ratio between the average coherent and 
incoherent signals can be represented as: 



(I c ) 1 A (n 2 +n) , 

{I ) Spin - T s) max(7,<5 p ) n l 

where 7 is the bandwidth of the final state. As long as the delay between the signal and the idler beams is smaller 
than 1/A, and in the absence of odd-order dispersion, we can assign P e sa 1. Taking into account that typically 
1 < gp(0) < 2, we see that the coherent signal is dominant not only at low photon fluxes (n -C 1, i.e. at the 
entangled-photons regime) but also at classically-high power levels, as long as both the pump and the final state of 
the interaction are narrower than the down-converted bandwidth: 

A » — t~tt <™> 

(P c ) max(7, dp) 

In the case of coincidence detection, the relatively long gating time T g of the electronic coincidence detection circuit 
makes the incoherent contribution to the coincidence counts rate Rc much larger: 

R*_ _J_ (n 2 +n) 

Rc lc 2TgA n 2 ' 1 ' 

Thus, for the coherent contribution to dominate in electronic coincidence-detection, one is restricted to very low 
photon fluxes (n <C 1). 



It is important to note that in this paper we took into account only two-photon interactions that result from cross- 
mixing of the signal and the idler fields and not from self-mixing of the signal with itself or the idler with itself. In 
both TPA and SFG, the cross-mixing term can be isolated spectrally if the signal and the idler are non-degenerate. 
In SFG, the cross-mixing term can also be isolated spatially if the down-conversion is non-collinear. However, in cases 
where the self-mixing term is indistinguishable from the cross-mixing terms (for example in the case of TPA with 
degenerate signal and idler fields, or if degenerate and collinear down conversion is considered) this has the effect of 
increasing the incoherent signal by a factor of two: 

I ic =*> 2P C (110) 

Naturally, the coherent signal is generated only by cross-mixing of the signal and the idler fields, and therefore is 
not affected by such self-mixing terms. 



Finally we note again that none of the effects described in this paper is directly related to squeezing. Even the non- 
classical linear intensity dependence is in fact independent of squeezing; since the coherent and incoherent signals are 
attenuated equally (quadratically) by such losses, this effect can be observed even in the presence of losses that would 
wipe out the squeezing properties completely Moreover, while the squeezing degree grows with n and is very small for 
n < 1, the linear term becomes less and less dominant as n grows, and is completely negligible at?i> 1. Furthermore, 
excluding the linear intensity dependence of the coherent signal, all the other effects considered in this paper are 
completely described within the classical framework. Indeed, such effects can be created by appropriately shaping 
classical pulses, so that they obtain similar anti-symmetric spectral phase correlations [52I. [53j. However, the precision 
of these correlations in broadband down-converted light can be many orders of magnitude higher than achievable by 
pulse-shaping techniques (42l . |43| . The unique properties of two-photon interactions induced by broadband down- 
converted light are therefore both interesting and applicable. 
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